Suppression of high pr hadrons in Pb+Pb collisions at LHC 
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Hadron production and their suppression in Pb+Pb collisions at LHC at a center-of-mass energy of 
^/sNN — 2.76 TeV are studied within a multiphase transport (AMPT) model whose initial conditions 
are obtained from the recently updated HIJING 2.0 model. The centrality dependence of charged 
hadron multiplicity dN c h/drj at midrapidity was found quite sensitive to the largely uncertain gluon 
shadowing parameter s g that determines the nuclear modification of the gluon distribution. We 
find final-state parton scatterings reduce considerably hadron yield at midrapidity and enforces a 
smaller gluon shadowing to be consistent with dN^/drj data at LHC. With such a constrained 
parton shadowing, charged hadron and neutral pion production over a wide transverse momenta 
range are investigated in AMPT. Relative to nucleon-nucleon collisions, the particle yield in central 
heavy ion collisions is suppressed due to parton energy loss. While the calculated magnitude and 
pattern of suppression is found consistent with that measured in Au+Au collisions at ^snn = 0.2 
TeV at RHIC, at the LHC energy the suppression is overpredicted which may imply the medium 
formed at LHC is less opaque than expected from simple RHIC extrapolations. Reduction of the 
QCD coupling constant a s by ~ 30% in the higher temperature plasma formed at LHC as compared 
to that at RHIC was found to reproduce the measured suppression at LHC. 
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Heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) 1-51 and recently at the Large hadron Col- 
lider (LHC) a, @] have revealed a new state of mat- 
ter comprising of strongly interacting quarks and glu- 
ons (sQGP) [7|, |8j. Primary evidence of this is provided 
by the observed suppression of high transverse momenta 
single hadron spectra [t| [l(| in central collisions relative 
to both peripheral and nucleon-nucleon collision. The 
suppression has been established as due to energy loss 
by the propagating partons in the plasma primarily by 
radiative gluon emission [ill E3 • Since the parton scat- 
terings occur at the early stage of the evolution in nuclear 
collisions, study of energy loss can probe the sQGP phase 
of the matter. In fact, the magnitude of energy loss is 
predicted to be strongly dependent on the parton density 
of the medium which reappears as soft hadrons [IH, [l3| . 

In addition to the final state parton energy loss, the 
jet quenching at moderate and high pr is also influenced 
by initial spatial distribution of partons, collective flow, 
and to the unknown nuclear shadowing of the parton dis- 
tribution. As the matter created in Pb+Pb collisions at 
\J s nn = 2.76 TeV is at about twice the density and 
probes parton distribution at a smaller momentum frac- 
tion x than at RHIC, analysis of the recent data for bulk 
hadron producti on jbH . ITa ] and high-pr hadron suppres- 
sion at LHC [lfjl . Il7| may provide crucial insight into the 
nuclear medium effects of parton shadowing and energy 
loss in the hot and dense QCD matter. 

While perturbative quantum chromodynamics 
(pQCD) can address only hard scatterings, formation of 
strongly coupled near perfect fluid as well as abundant 
soft particle production suggest a highly nonperturba- 
tive physics which is not yet well-established within 



QCD. Consequently models based on (non-)ideal hy- 
drodynamics [l8l - [20| . transport calculations I2l| - [23| . 
and transport/hydrodynamics hybrid models [24j have 
been developed. It was recently demonstrated in the 
HIJING 2.0 model [2|| that the larger uncertainties of 
the shadowing effects at RHIC [26] can be constrained 
from comparison of the measured charged particle 
density at midrapidity for the most central Pb+Pb 
collision at LHC. On the other hand, collision centrality 
dependence of bulk hadron observables should reflect 
the relative contribution to particle production from 
hard and soft processes. Thus a precise estimate of 
nuclear shadowing and detailed study of medium effects 
on particle production from soft to the hard scattering 
regime relies on systematic inclusion of various stages of 
dynamical evolution of matter. 

A MultiPhase Transport (AMPT) model [H| which 
combines the initial particle distribution from HIJING 
model [22} with subsequent parton-parton elastic scatter- 
ings via the ZPC parton cascade model and final hadron 
transport via ART allows a systematic study of hadron 
production and jet quenching. In this letter we shall 
investigate bulk charged particle production and jet sup- 
pression within the AMPT model modified to include the 
updated HIJING 2.0 version. In absence of control d+Pb 
data essential to calibrate the nuclear shadowing of ini- 
tial jet spectra, we shall use the centrality dependence 
of charged particle pseudorapidity density, dN^/dr], of 
the ALICE data in Pb+Pb collisions to provide a more 
stringent constraint on the gluon shadowing parameter 
s g which will be then employed to investigate jet sup- 
pression. 

In the two-component HIJING model [27j for hadron 
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production, nucleon-nucleon collision with transverse 
momentum p? transfer larger than a cut-off po leads to 
jet production calculable by collinearly factorized pQCD 
model. Soft interactions with pt < po is characterized by 
an effective cross section <r so ft- In the HIJING 2.0 model 
[25l | the Duke-Owens parametrization [28[ of the parton 
distribution functions has been updated with the modern 
Gluck-Reya-Vogt (GRV) parametrization (2j|. Since the 
gluon distribution at small momentum fraction x is much 
larger in GRV, instead of a fixed value for p — 2 GeV/c 
and (Tsoft = 57 mb (as used in HIJING 1.0), an energy 
dependent cut-off for po(y/s) and (J S oh(y/s) is used to fit 
experimental data on total and inelastic cross sections 
and hadron rapidity density in p + p/p collisions [25|. 

For the nuclear parton distribution function (PDF), 
HIJING employs the functional form 



f^x,Q') = AR*(x,Q*)f?(x,Q-') 



(1) 



where is the PDF in a nucleon. The nuclear modifi- 
cation factor of quarks and gluons (a = q, g) in HIJING 
2.0 parametrization are [25| 
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The impact parameter dependence of shadowing is taken 
as s a (b) — (5s a /3)(l — b 2 /R 2 ) that prohibits rapid rise 
of particle production with increasing centrality. Here 
Ra ~ A 1 / 3 is the nuclear size and s q = 0.1 is fixed 
by data from deep inelastic scatterings. From com- 
parison to the centrality dependence of charged parti- 
cle pseudorapidity density per participant pair of nu- 
cleons, (dN c h / dr/) / ((N pai t) /2) in Au+Au collisions at 
\J snn — 0.2 TeV, the gluon shadowing parameter in HI- 
JING 2.0 model has been constrained to s g — 0.17 — 0.22. 
Whereas a stronger constraint of s g = 0.20 — 0.23 has 
been obtained from the reproduction of dN^/drj AL- 
ICE data for the most central (head-on) Pb+Pb colli- 
sions at yj snn — 2.76 TeV. Albeit, HIJING ignores the 
final state interaction of particles, and such an estimate 
of s g is entirely from initial state effects. We shall show 
the influence of final state parton energy loss [Il|-[l3[ as 
well as hadronic rescatterings modify considerably the 
dN c h/dr) yield and thereby the magnitude of the initial 
state nuclear shadowing s g for gluon distribution. In the 
present study we shall use the string melting version of 
the AMPT where the hadrons from HIJING 2.0 are con- 
verted to their valence (anti)quarks and parton recom- 
bination is employed for hadronization. The coalescence 
of dominant soft partons and also relatively large num- 
ber of hard jets produced at LHC will thus contribute to 
the final charged hadron spectrum. In the Lund string 



fragmentation function f(z) oc z~ 1 (l — z) a exp(— brn^/z), 
where z is the light-cone momentum fraction of the gener- 
ated hadrons with transverse mass tut, we employ the de- 
fault HIJING values of a = 0.5 and b = 0.9 GeV~ 2 . Un- 
less otherwise mentioned, at both RHIC and LHC energy 
we consider the strong coupling constant a s — 0.33 and 
screening mass \i = 3.226 fm _1 *3(| that correspond to 
parton-parton elastic scattering cross section of a ~ 1.5 
mb in the parton cascade. 
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FIG. 1: Top panels: Pseudorapidity distribution for charged 
hadrons in Au+Au collision at RHIC energy of y / sjvjv = 0.2 
TeV and in Pb+Pb collision at LHC energy of ^/sjvjv = 2.76 
TeV for (0 - 5%) centrality. The AMPT model predictions 
are without any final state interactions as in HIJING (dashed 
line); with parton transport i.e. HIJING+ZPC (dashed- 
dotted line), and with further hadron transport as in AMPT 
(solid line). The solid circles are the measured values from 
the BRAHMS at RHIC [H and ALICE at LHC Q. Bot- 
tom panels: The rapidity distribution from AMPT for K + , 
7r + , p and net baryons, B — B, at the RHIC and LHC energies. 

Figure [1] shows the pseudorapidity distribution of 
charged hadrons in the 5% most central collision in the 
AMPT model in Au+Au at ^/s^ = 0.2 TeV and 
Pb+Pb at yj snn — 2.76 TeV. The results are with gluon 
shadowing parameter of s g — 0.15 (at RHIC) and s g — 
0.17 (at LHC) that are found to be in good agreement 
with the measured dN^/drj distribution from BRAHMS 
H [HI at yfsN~N = 0.2 TeV, and the dN ch /dr] (\r]\ < 
0.5) = 1601 ± 60 from ALICE QJ] at v^vw = 2.76 TeV. 
In absence of final state partonic and hadronic scatter- 
ings, which is basically the HIJING 2.0 model predicts 
dN ch /dr] (\r)\ < 0.5) = 706 ± 5 and 1775 ± 3 at RHIC 
and LHC, respectively. In subsequent parton cascade 
(i.e. HIJING plus ZPC), energy dissipation and redistri- 
bution into the transverse flow via partonic scatterings 
lead to a reduction of charged particle multiplicity by 
surprisingly a similar amount of ~ 15% at both RHIC 
and LHC. Though the partonic density at LHC is about 
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twice than at RHIC, this nearly equal suppression of yield 
after parton cascade reflects the interplay between hard 
and soft processes via a delicate balance between col- 
lective flow, gluon shadowing and jet multiplicity all of 
these are larger at LHC than at RHIC. Finally, subse- 
quent hadronic scatterings (dubbed as AMPT) from the 
less dense phase leads to a smaller decrease of particle 
multiplicity. Fig. Q] further shows that final state scat- 
terings essentially smoothen out the dip at 77 = (due to 
Jacobian) in HIJING to a nearly flat dN^/drj distribu- 
tion around mid-rapidity. Such a weak pseudorapidity- 
dependence in dN r h / dy at 77 < 2 has also been observed 
in both the BRAHMS [3l| and CMS data 

The rapidity distribution of pions, kaons, antiprotons 
and net baryons are displayed in Fig. Q] at ^/snn = 0.2 
and 2.76 TeV. With more than an order of magnitude in- 
crease in energy at LHC, the rapidity distribution of the 
produced hadrons becomes wider by ~ 55% and thereby 
dN c h/drj at midrapidity increases by ~ 2.4 compared to 
the top RHIC energy. While the net-baryon density is 
found to decrease by ~ 35% from ^/snn = 0.2 TeV to 
2.76 TeV, the antibaryon to baryon ratio at these RHIC 
(LHC) energies are found to be p/p = 0.71(0.88), A/A = 
0.75(0.95), 3/E = 0.83(0.99) and CI /SI = 0.89(1.00). The 
yield ratios from the AMPT at RHIC are consistent with 
the feed down corrected measured values [l|, ||J within 
the systematic errors. Enhanced meson production and 
slight decrease in the strangeness density at LHC result 
in the ratios at midrapidity of p/ir + = 0.091 (0.088) and 
K+/tt+ = 0.17 (0.15) at the RHIC (LHC) energies con- 
sidered here. 

In Fig. [3] we present the charged particle pseudorapid- 
ity density per participant pair, (dN c h/ drf) / \(N par t) /2), 
as a function of centrality of collision characterized by 
average number of participating nucleons (V part ). The 
AMPT calculation are for Au+Au collisions at ^/snn — 
0.2 TeV with a range of gluon shadowing parameter 
s g = 0.10 - 0.17 and for Pb+Pb at y^T^F = 2.76 TeV 
with s g = 0.16 — 0.17. With this choice of the gluon shad- 
owing parameter, the centrality dependence of charged 
particle multiplicity agrees well within the experimental 
uncertainty seen in the BRAHMS [3ljj and PHENIX [H[ 
data at RHIC. Due to abundant jet and minijet produc- 
tion at LHC, the ALICE multiplicity data for Pb+Pb 
collision is quite sensitive to nuclear distortions at small 
x and provides a much stringent constraint on the gluon 
shadowing of s g ~ 0.17. It may be mentioned that the 
estimated values of s g in AMPT are consistently smaller 
than in HIJING 2.0 model [25| which underscores the im- 
portance of final state interactions in precise estimation 
of the nuclear shadowing of partons that in turn should 
also influence the hard observables. 

The study of bulk hadron production when coupled 
with that for hadron spectra provide crucial information 
of the parton- medium interactions where high-py partons 
suffer energy loss that are transported to produce soft 
hadrons. To quantify such a suppression of hadrons at 
high pt due to medium effects in heavy ion collisions, the 
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FIG. 2: Charged hadron multiplicity density dN c h/dri at mid- 
rapidity per participant nucleon pair as a function of average 
number of participants {N pSlI t). The results are from AMPT 
calculations (triangles) obtained with gluon shadowing pa- 
rameter s g = 0.10 — 0.17 in Au+Au collision at y^ijvjv = 0.2 
TeV (top panel) and with s g — 0.16 — 0.17 in Pb+Pb collision 
at t/snn = 2.76 TeV (bottom panel) as compared with the 
data (circles) from BRAHMS and PHENIX Q at RHIC 
and ALICE [3 at LHC. 



nuclear modification factor 



Raa(pt) 



d N / dt] dpT 
(N colx ) d 2 NPP/d V d Pl 



(3) 



is used which is the ratio of particle yield in heavy ions 
(A + A) to that in p + p reference spectra, scaled by 
the total number of binary nucleon-nucleon (NN) col- 
lisions (iVoou) = (Taa) Cfael- I n absence of initial and 
final state nuclear medium effects Raa(pt) — 1 by con- 
struction. The nuclear thickness function (Taa) and the 
inelastic NN cross section a^ei are calculated within the 
HIJING 2.0 model that uses Glauber Monte Carlo simu- 
lation for distribution of initial nucleons with a Woods- 
Saxon nuclear density. The energy dependent soft in- 
teraction cross section a so [t(\/s) in HIJING 2.0 enforces 
o'iniT t0 be about 42 and 64 mb at ^/sWlv = 0.2 and 2.76 
TeV, respectively. However, at low pt regime dominated 
by soft particle production, the scaling by the number 
of nucleons suffering at least one inelastic collision, i.e 
Vp ar t, is more appropriate. 

Figure [3] shows the inclusive charged hadron pt spec- 
tra at midrapidity in the AMPT for p + p collisions and 
for central (0 — 5%) and peripheral Au+Au collisions at 
^Jsnn — 200 GeV (left panel) and in Pb+Pb collisions at 
\J s nn = 2760 GeV (right panel). The results are for ini- 
tial parton distribution with gluon shadowing s g = 0.15 
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FIG. 3: Invariant hadron production spectrum in p + p and 
Au+Au collision at ^snn = 0.2 TeV (left panel) and in p+p 
and Pb+Pb collision at s/snn = 2.76 TeV (right panel). The 
results are from AMPT calculations in p + p (dashed lines) 
and heavy ion (A + A) (open symbols) collisions with gluon 
shadowing parameter s g = 0.15 (0.17) at RHIC (LHC). The 
measured spectrum are for Au+Au (solid symbols) and p + p 
non-single-diffractive interaction (star) by STAR [9( at RHIC 
and for Pb+Pb (solid symbols) by ALICE Q]| at LHC. The 
p + p reference spectrum at y / ijvjv = 2.76 TeV (solid lines) 
is the ALICE interpolation normalized to LO pQCD which 
is shown as scaled by average number of binary collisions, 
(iVcoii), corresponding to the centrality classes. 



(0.17) at RHIC (LHC) energies that have been fixed from 
the centrality dependence of N c h data. In p+p collisions, 
the pr spectra from the model exhibit the LO pQCD 
based power law behavior at px > 5 GeV/c which is 
in overall good agreement with the STAR data [9(. At 
y?SNN = 2.76 TeV, we however find the calculated yield 
from p + p overpredicts at p T > 6 GeV/c that obtained 
by ALICE [l6j from interpolation of pp spectrum mea- 
surements at yj sjv n = 0.9 and 7 TeV to estimate the 
suppression Raa- For peripheral heavy ion collisions the 
AMPT spectra are consistent with that measured at both 
RHIC and LHC energies. On the other hand, the pr 
distributions for central collision show marked deviation 
from power law function and are clearly suppressed es- 
pecially at moderate px = 4 — 11 GeV/c due to medium 
modification. Though the AMPT spectra from central 
collisions describes the RHIC data quite well, it is how- 
ever much softer than the ALICE data at pr > 2 GeV/c. 
This possibly stems from enhanced energy loss of partons 
in a much denser medium that is generated from melting 
of strings to their valence quarks and antiquarks in the 
QGP medium. 

The nuclear modification factor Raa for charged 
hadrons is shown in Fig. [4] for central and peripheral 
Au+Au collision at RHIC (top panel) and in Pb+Pb 
collisions at LHC (bottom panel). For central collisions 
at both energies, Raa(pt) is less than unity which im- 
plies appreciable suppression of charged hadrons rela- 
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FIG. 4: Nuclear modification factor Raa for charged hadrons 
and neutral pions as a function of pr in central and peripheral 
Au+Au collisions at y / sjviv = 0.2 TeV (top panel) and Pb+Pb 
collisions at ^/sWn = 2.76 TeV (bottom panel). The AMPT 
model predictions are compared to the data from STAR 0] 
and PHENIX [hJ at RHIC and from ALICE Q]| at LHC. The 
AMPT results are with strong coupling constant a s = 0.33 
at RHIC and LHC and with a s = 0.24 for central collisions 
at LHC. The histograms is the systematic error band due to 
different interpolation procedure used in earlier estimates by 
ALICE for the baseline p + p spectra. 



tive to NN reference. The model calculations, with nu- 
clear shadowing parameter s g — 0.15 constrained from 
dNch/dr) data in Au+Au collisions, describes the magni- 
tude and pattern of the RHIC suppression data @. It 
is seen that Raa increases gradually with px reaches a 
maximum of Raa — 0.7 at pr — 1.8 GeV/c, then it de- 
creases with further increase of pr and saturates there- 
after to about 0.2 at pr > 7 GeV/c. The success of 
AMPT at y/sWN = 0.2 TeV thus suggests that the ini- 
tial state shadowing of pQCD jet spectra, the final state 
scattering and the parton energy loss is consistent with 
the formation and evolution of the medium at RHIC en- 
ergy. 

At 70 — 80% centrality Pb+Pb collisions at ^/smn — 
2.76 TeV, the Raa for charged hadrons is nearly con- 
stant at about 0.7 over a large pr range as seen in both 
the ALICE data [3 and AMPT model calculations. At 
this peripheral collision, the QGP even if formed, should 
have a small volume and short lifetime. In central Pb+Pb 
collisions at LHC, the rise and fall pattern exhibited by 
Raa up to pr ~ 6 GeV/c is similar to RHIC. However, 
as evident from ALICE data, the suppression of charged 
hadrons at low px is somewhat larger, and Raa reaches 
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a minimum of 0.14 around 6-7 GeV/c. The earlier es- 
timates with large errors as shown by histogram is due 
to interpolation procedure used by ALICE to obtain the 
baseline p + p spectrum. With the recently measured 
spectrum in p + p collision at ^/snn — 2.76 TeV [Bj], 
the measured Raa drops but remains well within the 
systematic error bands which is also consistent with the 
CMS data [13 ■ In contrast to ALICE data, the AMPT 
calculations show even more pronounced suppression at 
Pt > 2 GeV/c due to significant quenching of the hard- 
scattered partons. Within the coalescence mechanism 
for hadronization in AMPT, though the peak positions 
and the subsequent decreasing pattern of Raa are simi- 
lar to the measured RHIC and LHC data, the minimum 
is found to be at 0.09 at px ~ 6 GeV/c. The subse- 
quent rise of Raa (compared to nearly constant value at 
RHIC) essentially stems from harder unquenched pQCD 
jet spectra at LHC and found to have similar slope as in 
the data. 

In Fig. [4] we also show the Raa for neutral pions for 
central collisions. As seen in charged hadrons, the Raa 
for 7T° exhibit a similar but a gradual rise and fall pattern 
at intermediate px (1.8 < pt < 5 GeV/c) and thereafter 
saturates (rises) with increasing pt at RHIC (LHC) ener- 
gies. Both the calculation and PHENIX data show that 
relative to charged hadrons, the 7r°s are more suppressed 
by as much as ~ 45% at the intermediate px- How- 
ever, at pt ^ 5 GeV/c the magnitude of suppression are 
nearly same for neutral pion and charged hadrons. The 
larger Raa for charged hadrons compared to neutral pi- 
ons can be explained as due to large baryonic (protons 
and antiprotons) yield produced from parton coalescence 
used for hadronization [33|, [3J| . In fact we find the in- 
variant yield of pions and protons become comparable at 
p T 2 — 4 GeV/c. At pt ^ 6 GeV/c, as pions are the 
most abundant particles and moreover the parton spectra 
become gradually flatter with increasing px, coalescence 
of hard partons is seen in AMPT to predict in nearly 
identical suppression Raa for pions and hadrons. 

The significant suppression in AMPT much below than 
the ALICE data suggests that the medium with more 
than a factor of two larger parton density than RHIC is 
in fact more transparent than expected. Attempt to in- 
crease Raa at high px by decreasing the shadowing s g 
only result in an enhanced bulk (soft) hadron produc- 
tion and thus disagree with the centrality dependence 
of dN ch /dii data shown in Fig. [2] In fact, the WDGH 
jet energy loss model [35j that has been constrained to 
fit the RHIC suppression data severely underpredicts the 
central Raa value of ALICE. 

On the other hand, we note that the above suppres- 
sion was calculated in the AMPT with same values of 
QCD coupling constant a s — 0.33 and screening mass 
fi = 3.226 fm- 1 at both RHIC and LHC. Perturba- 
tively, the screening mass depends on temperature as 
/i = gT with g = y/4ira s [36j . The parton-parton elas- 
tic scattering cross section used in AMPT then reduces 
to a 9na 2 /(2fi 2 ) 9a s /(8T 2 ). Since hydrodynamic 



model analysis of RHIC/LHC data indicate [33] only 
about 10% viscous entropy production, the initial en- 
tropy density can be approximated to final particle 
multiplicity by the scaling 38] 
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where A± is the transverse area of the collision zone. 
The proportionality constant for entropy rapidity den- 
sity, dS/dy, to dN c h/dy conversion was taken from Refs. 
38 40]. For a QGP characterized by massless gas of 
light quarks and antiquarks, Si ~ 4ej/(3Tj) with en- 
ergy density (21/30)7r 2 7; 4 . This allows to esti- 
mate the initial temperature Tj and thereby the par- 
ton scattering cross section a from the measured par- 
ticle yield. For the 5% most central Au+Au and Pb+Pb 
collisions at ^/snn = 200 and 2076 GeV, the measured 
dN ch /dy w 687 [32| and 1601 Q result in T t w 320 
and 436 MeV respectively, at a proper time r, = 1 
fm/c. With the above choice of a s = 0.33, the esti- 
mated cr « 9a s /(8T 2 ) « 1.4 mb at RHIC is incidentally 
close to the value employed in AMPT that reproduces 
the RHIC suppression data shown in Fig. S) In con- 
trast, the higher temperature T{ at LHC enforces a much 
smaller a ~ 0.76 mb. Alternatively, if the screening mass 
remains constant at (i = 3.226 fm _1 from RHIC to LHC, 
such a small a is then consistent with a s ~ 0.24 at LHC. 
With this reduced a s , we show in Fig. 2] (open circles) 
the AMPT results of Raa of charged hadrons in cen- 
tral Pb+Pb collisions at ^/sWn = 2.76 TeV. The good 
agreement with the ALICE suppression data is a clear 
indication of thermal suppression of the QCD coupling 
constant due to higher temperature at LHC compared to 
that at RHIC. It may however be mentioned that instead 
of an average value, the strong coupling could have a 
temperature dependence a s (T) during the plasma, evolu- 
tion [4l|. Further, the AMPT model calculations invoke 
purely elastic collisional energy loss, the effects of inelas- 
tic scatte rings v ia medium-induced radiative parton en- 
ergy loss [ll|, [l2|, l23j, l35J could still pose a serious theoret- 
ical challenge to understand the underlying energy loss 
mechanism especially at the LHC energy regime. 

In summary, we study the nuclear medium effects on 
hadron production over a wide range of px in Pb+Pb 
collisions at the LHC energy ^/snn — 2076 GeV. For 
this purpose we use the AMPT model which is updated 
to include the HIJING 2.0 version for initial conditions 
for parton distribution. We find final-state parton scat- 
terings reduce significantly the hadron multiplicity at 
midrapidity that enforces smaller gluon shadowing for 
agreement with the ALICE data for charge particle yield 
at various centralities. With such a constrained parton 
shadowing, we find that parton energy loss in AMPT 
describes quite well the magnitude and suppression pat- 
tern of hadrons in both central and peripheral Au+Au 
collisions at the RHIC energy ^/s NN — 200 GeV. With 
the same strong coupling constant a s = 0.33, the model 
however predicts larger jet quenching relative to ALICE 
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data for central Pb+Pb collisions at LHC. A reduction of 
a s by ~ 30% in the higher temperature plasma formed 
at LHC was found to describe the measured suppression. 
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